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I. INTRODUCTION 

The dramatic progress in laser technology has enabled 
a novel area of studies exploring laser-matter interac- 
tions at ultra- high intensity pj. The intensity level of 
2 x 10 22 W/cm 2 has recently been achieved and two 
projects [iilH aiming at intensity levels up to 10 26 W/cm 2 
have been supported and are under way. Furthermore, 
several original proposals have been suggested e.g. 0- 
0], which reach even higher intensities with almost the 
present level of technology. One of the key phenomena 
of laser-matter interactions, that probably dominates at 
ultra-high intensities of our interest, is the occurrence 
of QED cascades [3, l^-Hoj . These cascades (also called 
avalanches, or showers) are caused by successive events of 
hard photon emissions and electron-positron pair photo- 
production by hard photons. As predicted in [lfl | based 
on qualitative estimations, the cascades may arise as soon 
as the laser field strength exceeds the threshold value of 
E+ = aEs, where a = e 2 /hc « 1/137 is the fine struc- 
ture constant and E$ = m 2 c 3 /eh = 1.32 x 10 16 V/cm is 
the characteristic QED field. Such a field strength corre- 
sponds to an intensity of ~ 10 25 W/cm 2 . 

Previously QED cascades have been observed and 
studied as a part of Extensive Air Showers (EAS) in the 
context of the passage of ultra-high energy particles, that 
orig inate from Cosmic Rays, through the atmosphere 
[Ill4l3| . However, similar processes can be observed in 
the external electromagnetic field as well. In this case, 
Brcmsstrahlung is replaced by the non-linear Compton 
scattering and Bethc-Hcitlcr process is replaced by the 
non-linear Breit- Wheeler process. The latter processes 
are well studied both theoretically (l3 - [l9| and in laser 
experiments [20| and are pro bably of great importance 
for astrophysics (see, e.g., |21|). 

An important novel distinctive feature of the cascades 
in the ultra-strong laser field, compared to situations ever 
studied previously, is that the laser field is not only able 
to be a target for ultra-relativistic electrons and hard 
photons, but can also accelerate the charged particles to 



ultra-relativistic energies. 

As a result, the cascades can be produced even by 
initially slow electrons or positrons, if they were some- 
how injected into the strong field region. Moreover, the 
mean energy of the particles is no longer decreasing in 
the course of the cascade development due to its redis- 
tribution among the permanently growing number of the 
created particles, but rather is restoring at the expense 
of the energy extracted from the laser field. This must 
lead to a vast increase of the cascade yield, as compared 
to the cascades in media or in strong magnetic fields. In 
this case the cascade multiplicity would be restricted ei- 
ther by the duration of stay of the particles in the focal 
region of the laser field, or even, under more extreme con- 
ditions, by the total energy stored in the laser field. In 
the latter case the focused laser pulses would be depleted 
by cascade production. 

As it will be explained in more details below, the 
restoration mechanism works if the particles can be ac- 
celerated transversely to the field. It was conjectured 
H |T(| on the basis of qualitative analysis for the model 
of a uniformly rotating electric fields, that this may be 
indeed the case. 

In the EAS theory, the ID approximation is often used 
because spreading in transverse direction is inessential for 
ultrarelativistic particles and has no significance for that 
problem. Besides, the cascade equations can be solved in 
this case analytically within the ultra-relativisti c ap prox- 
imation by means of the Mellin transform [l2|, [l3|]. The 
results of such analytic theory are in good agreement with 
both experiments [111 [l2j and direct Monte- Carlo simu- 
lations [22[ • The attempts to treat the cascades in strong 
magnetic fields on similar grounds are also known [23| . 
However, though the ID approximation remains valid, 
the cascade equations can not be simplified via the Mellin 
transform unless some further approximation is made. 
According to Monte-Carlo simulations [24| , the resulting 
analytical approach works much worse here than in the 
case of ID approximation for EAS. In our case of cas- 
cades arising in a laser field, the structure of the cascade 
equations (see Appendix [X| is the same as for the mag- 
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nctic field, but it is impossible to incorporate restoration 
mechanism within the ID approximation in momentum 
space. This means that our problem is essentially two- 
or three dimensional. 

In this work we report on the first results of the Monte- 
Carlo simulations of cascades produced by initially slow 
electrons in a uniformly rotating homogeneous electric 
field. Such a field can be obtained practically at the 
antinodes of a standing electromagnetic wave. The choice 
of the field model is uniquely specified by the existence of 
reasonable qualitative estimations for scaling of the basic 
cascade characteristics for this particular case Our 
goal was to prove explicitly the existence of the restora- 
tion mechanism and to test the estimations [l(| by direct 
numerical simulations. 

The paper is organized as follows. In Sec. [HI which can 
be considered as a technical introduction, we review and 
collect the known information on the elementary quan- 
tum processes: single photon emission by electrons and 
pair creation by hard photons in strong fields of arbitrary 
configuration. Though this information is not completely 
new, it is of essential importance for our presentation and 
is spread among the literature on the subject. After that, 
in Sec. Mil we present the reasoning in favor of the en- 
ergy restoration mechanism for cascades in electromag- 
netic fields. In Sec. IIVI we formulate the assumptions of 
our model, present the details of our Monte-Carlo rou- 
tine and discuss the results obtained by numerical sim- 
ulations. These results are compared to the known es- 
timations. Summary and discussion is given in Sec. [V] 
Finally, in Appendix [XJ we discuss the cascade equa- 
tions for our problem and explicitly demonstrate that, 
contrary to the recent doubts [25j], the approach we use 
takes proper account of radiation friction by ultrarela- 
tivistic electrons. 



II. QUANTUM PROCESSES WITH 
HIGH-ENERGY PARTICLES IN A STRONG 
ELECTROMAGNETIC FIELD 

General properties of radiation of ultrarclativistic par- 
ticles are well known [2(| . Due to the relativistic aberra- 
tion effect the momenta of the products of any decay of 
an ultrarelativistic particle are directed within the nar- 
row angle A9 ~ 7 _1 with its momentum, where 7 = 
y/l + (p/mc) 2 . Thus, radiation of a charged ultrarela- 
tivistic particle is visible at a point of observation only 
for a short period of time r during which its momentum 
turns through the angle of the order A8. The momen- 
tum turning angle can be estimated by A8 ~ eF±r / mcy , 
where F± denotes the characteristic value of the trans- 
verse component of the field. Thus, r ~ mc/eF±. The 
characteristic frequency f2 of classical radiation can be 
most simply estimated in the proper reference frame of 
the particle, by transition to which the duration r trans- 
forms into t' = r/7. Thus, ft' ~ r' 1 ~ eF±j/mc. In 
the laboratory frame, due to the Doppler up-shift, we 



would thus have f2 ~ (eF±/mc)-f 2 . Such a scaling with 7 
is typical for congeneric problems and arises e.g. in the 
theory of synchrotron radiation. 

The parameter x = hVl/^mc 2 ) = F±^/Es (27[, being 
the ratio of the classically estimated mean energy of an 
emitted photon to the energy of the radiating particle, 
determines whether the process of radiation is controlled 
by classical electrodynamics or QED. Namely, if \ -C 1, 
then quantum recoil is inessential and radiation is classi- 
cal, whereas if % > 1 then it must be quantum. The pa- 
rameter \ is Lorentz and gauge invariant and is precisely 
defined as x — eh/(m 3 c 4 )^y— (F^p") 2 , where F^ v is the 
strength tensor of the field and p^ 1 is the 4-momentum 
of a particle. In what follows, we assume that F <C Es- 
On the other hand, we assume that the field is of rela- 
tivistic strength in the sense that the dimensionlcss field 
amplitude ciq = eW—A^A^ / (mc) ^> 1, where A^ is the 
4-vector of the field potential. The latter means, in par- 
ticular, that it varies on the scale that exceeds r and thus 
can be considered constant with respect to the decay pro- 
cesses. 

Two different theoretical approaches have been devel- 
oped in order to study photon emission by ultrarelativis- 
tic (7 3> 1) charged particles in electromagnetic fields 
of ultrarelativistic (ao 3> 1, X ~ 1) but still subcritical 
(F -^i Es) intensities. Nikishov and Ritus (NR) have cal- 
culated the appropriate quantum amplitudes in terms of 
Volkov solutions in a constant crossed field (E 2 — H 2 = 
and E ■ H = 0) of arbitrary strength [27|. Their re- 
sults can be applied directly to our problem because, as 
they have pointed out especially, under the abovemen- 
tioned conditions any field looks locally as constant and 
crossed, the latter in the sense that both field invariants 
(E 2 - H 2 )/E 2 and E ■ H/E 2 are much less thanv 2 . 
The other approach by Baier and Katkov (BK) fl9H27| is 
based on the observation that the motion of a particle in 
between two acts of photon emission (which corresponds 
to free lines of Feynman diagrams) can be considered 
classically if F -C Es, so that all the relevant quantum 
corrections are reduced to quantum recoil and, possibly, 
to field-spin interactions. In a sense, the BK approach 
is equivalent to the replacement of the aforementioned 
Volkov solutions by the localized wave packets moving 
along the classically prescribed trajectories. Both ap- 
proaches are based essentially on the same approxima- 
tions and provide the same results for the energy spectra 
of emitted photons. 

The energy distribution of the probability rate for pho- 
ton emission by ultrarelativistic electrons in an electro- 
magnetic field is given by [ljl Qjl [27j 
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where x = (x 7 /XeX'e) 2/3 , Ai(x) = (1/tt) J °° cos(£ 3 /3 + 
£x) d£ is the Airy function, e 7 and e e are the energies of 
the emitted photon and the initial electron, respectively. 
Xe, x' e = Xe ~ X 7 an d X 7 (0 < Xi < Xe) are the di- 
mensionless quantum parameters for the electron before 
and after emission, and for the emitted photon, respec- 
tively. In terms of the field strengths, this parameter is 
represented as 
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where e and p are the energy and the momentum of the 
corresponding particle. 

Note that expression ([1]) suffers from the infrared sin- 
gularity at e 7 — > 0. However, in our constant field ap- 
proximation this singularity dW ra d{e-y)/de-y = 0(e 7 2 ^ 3 ) 
is weaker than the usual O(eZ^) scaling of the infrared 
behavior of perturbative QED [23, H^|. In particular, the 
total radiation probability rate is infrared convergent in 
our approximation. The infrared sector, however, is not 
important for the parameters considered in our paper, 
because most of the emitted radiation are found to have 
much larger frequencies than the frequency of the driving 
field. 

The energy distribution of the probability rate for di- 
rect pair creation by hard photons (e 7 S> mc 2 ) is given 
by EEHg 
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where the indices "7" and "e" this time refer to the initial 
photon and to the created electron, respectively. For 
the created positron, we have x'e ~ X-y ~ Xe (0 < x 7 < 
Xe). Formula (J3]) is completely symmetric with respect 
to electron and positron remaining unchanged under the 
replacement Xe x' e - Similarity between formulas ([1]) 
and J3|) is explained by the fact that these two processes 
arc related by the cross-symmetry [27j |. 

The total probability rates for both processes 
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cannot be written in terms of known special functions and 
should be obtained by numerical integrations. However, 
they allow simple asymptotic expressions in the limits of 
small and large Xe, X 7 ? respectively. Namely, we have 



W rad « 1.44— x e , Xe < 1, 



(5a) 




FIG. 1: The sign of Xe(i) along the particle trajectory at 
t — to in different zones of the poxPoy plane. The shaded zones 
correspond to acceleration (increase of Xe) of positrons and 
deceleration (decrease of Xe) of electrons. The non-shaded 
zones correspond to the vice-versa situation. 
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Eqs. (|5a|) . (|6a|) in these formulas describe the quasiclas- 
sical regime. For small values of the quantum parameter 
X 7 , the probability rate for pair photoproduction W cr is 
suppressed exponentially, in accordance with the essen- 
tially quantum nature of this process. At the same time, 
W rad remains 0(h^ 1 ) 1 thus providing a finite classical 
limit for the mean radiated intensity I ra d- As for the 
limit of large x, both rates (|5b[) and ()6b[) differ only by 
a numerical factor of the order of unity. 

Given the energy and the momentum of the electron 
before emission, Eq. ([T]) determines the probability dis- 
tribution for the energy e 7 of the emitted photon. Under 
our assumption 7 e ^> 1 , the momentum of this photon is 
given by p 7 = (e 7 /p e )Pe- The energy and the momen- 
tum of the electron after emission should be determined 
from the conservation laws. In the electromagnetic back- 
ground, they are of the form p£ + q^ 1 = p' e + pf^, where 
qf 1 is the four-momentum extracted from the field. The 
exact value of this g M essentially depends on the global 
structure of the field. This is because the whole space- 
time contributes to the integrals in the QED matrix ele- 
ment that yield delta-functions expressing the conserva- 
tion laws. For example, in a crossed constant field with 
the Poynting vector directed along the z-axis the con- 
served quantities arc s/c — p Zl p x and p y . In a constant 
electric field, the canonical momentum is conserved. In 
a constant magnetic field, directed along the z axis and 



4 



Non-perturbed 




Initial position 



-1 0123456 

x, c/uj 

FIG. 2: Spatial picture of the formation of the cascade ini- 
tiated by a positron in the homogeneous uniformly rotat- 
ing electric field (obtained by a Monte-Carlo simulation with 
ao = 2 x 10 3 and ftw = leV). Legend: Trajectories of elec- 
trons and positrons are shown as black and gray curves, re- 
spectively. The hard photons which have created pairs during 
the simulation time are shown as the dashed lines. The tra- 
jectory of the primary positron ignoring any QED processes 
is plotted as the thick light gray curve. 

with symmetric gauge, the conserved quantities are the 
number of the Landau level, the angular momentum, and 
p z . Nevertheless, for ultrarelativistic particles there is ac- 
tually no difference between these possibilities, and cither 
of them can be adopted with the accuracy of our approx- 
imation. The reason is that q < eFr < mc <C p e ,p' e ,p^. 
In particular, we can assume p' e = p e ~p 1 . The same ar- 
gument can be applied to pair creation by hard photons 
as well. 

In addition to one-photon emission and direct pair 
photoproduction reviewed above, there exist more com- 
plicated higher-order processes, such as e.g. the two- 
photon emission e~ — > e~77 or the trident process 
er — > e~e~e + . Their specific feature is that the interme- 
diate particle is off the mass shell, i.e. is virtual. How- 
ever, in strong field situations of our interest the two-step 
processes dominate [1, [^, [2(| ■ For this reason, we do not 
consider higher-order processes in the sequel. 



III. BASIC ESTIMATIONS FOR CASCADE 
PRODUCTION IN A ROTATING ELECTRIC 
FIELD 

Since there is no difference whether an electron or a 
positron initiates a cascade, we assume in this section 



that our cascade is initiated by a positron (e > 0). 

Consider a positron in a homogeneous, uniformly ro- 
tating electric field 

E(t) = {E cos ujt,E smut}. (7) 

The equation of motion 

p(t) = eE(t), (8) 

with the initial condition p(to) = pn, can be easily solved: 

p x {t) = Pox + mcan(smujt — sinwto), 

p y (t) = p y o — mcao{cos uit — cosuito)- (9) 

Here, ao = eEn/mujc is the dimensionless field amplitude. 

Let us assume first that the positron is at rest (po = 0) 
initially (to = 0). Equations ([7]) and ([9]) show that the 
energy and the quantum parameter x of the positron for 
this case depend on time as 



e e (f) = mc 2 \j 1 + 4a 2 , sin 2 —, (10a) 
, . ehEn / „ , cut , , , 

Xe« = ^V 1 + 4a ° Sin V (10b) 

Both quantities are increasing initially. They are os- 
cillating with the period 2tt/lj of the rotation of the 
field. The amplitudes of these oscillations, e m « 2mc 2 ao, 
Xm ~ 2ao(£'o/ Es) = 2{hu> / mc 2 )a^ are proportional to 
ao and a 2 ,, respectively, and are quite large under our 
basic assumptions. For example, Xm approaches unity 
already at ao ~ ao c = 500 for an optical rotation fre- 
quency of fko = leV. This corresponds to the field 
strength Eq ~ IQ~ 3 E S ~ 10 13 V/cm and the intensity 
10 24 W/cm 2 . Since Xm ~ 1 under such conditions, the 
positron, according to the preceding section, is able to 
emit a hard photon with Xi ~ Xe ~ lj which, in turn, 
can create an electron-positron pair. However, at such in- 
tensities a new generation of pairs is typically produced 
on the time scale 7r/w, and the whole pair generation pro- 
cess may be rather sensitive to peculiarities of the field 
model. As we discuss below, stable cascade formation is 
expected at higher intensity levels. 

The formulas (|10aP and (|10b[) become especially simple 
for stronger fields ao S> ao c , because in this case the value 
Xe ~ 1 is being reached within just a small fraction t acc 
of the rotation period. Namely, we have 

eJt) w eE ct, — <t<i, (11a) 
wao w 

1 / En \ 2 mc 2 UJ o 1 1 / , s 

2 \Es J n tjy/oo ijj 

Eq. (|lla|) is easy to understand, because initially the 
positron is accelerating almost along the field. In or- 
der to understand Eq. (|llb[) better, let us note that in 
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FIG. 3: Comparison of the cascade profile obtained with our 
code [black thin line], with a code applying an alternative 
event generator [circles] and from previous independent sim- 
ulations [thick gray line, see Fig. 5 in H3]]- Depicted is the 
number of pairs with energy e > 10 _3 £o versus the elapsed 
time. The simulation parameters are £o = lOOGeV and 
E /Es = 0.2. 



the case p(0) = 0, according to Eqs. Q and ©, the 
momentum of the positron constitutes exactly the angle 
Lot/ 2 with both the instant direction of the field and the 
x-axis. Intuitively, this is because due to its inertia the 
particle does not follow the rotation of the field precisely. 
As a consequence, the transverse component of the field 
with respect to momentum of the particle increases as 
E± = E sm(ujt/2) « E ujt/2. Since xS) ~ E ±-l/ E S, 
we arrive immediately at Eq. (|llb[) . Qualitatively, the 
same growth of the energy and the parameter x with 
time has been observed for generic field configurations 

As it follows from Eq. (|llb[) , the quantum parameter 
Xe becomes of the order of unity over the period of time 



Eq. © and consider the general initial condition. It is 
easy to see that in the case H = the sign of the deriva- 
tive of the quantity @ is determined completely by the 
expression —e(p-E){jp-E). The zones in the po-plane in 
Fig. [TJ where Xe > is valid for positrons at time t = to 
are shaded. In the shaded areas x& can be expected to 
increase for some time. The time io can be identified 
with the creation time of a new pair. 

Since the momentum of a primary positron is confined 
to the shaded region to the right, and the new secondary 
particles are created with momenta parallel to the mo- 
mentum of parental particle in our approximation, we 
see that momenta of the secondary particles also lie in 
the shaded sector. Thus, the newly created positrons are 
accelerating with x(t) > 0, while the newly created elec- 
trons are initially decelerating (x < 0). However, they 
are quickly turned back by the field and also get acceler- 
ated. The same will be true for successive pair creation 
processes as well. 

In spite of the complexity of the picture of the cas- 
cade development (see Fig. [2]), some general estimations 
for it can nevertheless be obtained [Tq] in the high-field 
limit /i 1. The idea is that, due to similarity of the 
rates (|5b[) and (|6b[) . and also because the variation of the 
angles between the momenta of all particles (positrons, 
electrons, and photons) and the field is determined by 
the same temporal scale there is actually no need 
to distinguish between all three sorts of particles. So, an 
order of magnitude estimation can be provided within 
the model of a simple doubling chain process. 

Let us denote by t e the typical lifetime for electrons 
and positrons with respect to hard photon emission. The 
same quantity up to an order of magnitude defines the 
lifetime of photons with respect to pair creation. The 
lifetime t e , together with the typical energy and the value 
of quantum parameters of the particles, as well as with 
the angle between their momenta and the field, can be 
estimated as [Toj 



t„ 



(12) 



Here we have introduced a new dimensionless field inten- 
sity parameter fx = E/E*, E* = aE$ « Sg/137, which is 
appropriate for the cascade problem [l(| . The parameter 
fx is related to the commonly accepted parameter ao by 
fx = (hu/amc 2 )ao- According to Ref. [1J|, the cascades 
can be caused by initially slow particles if fx > 1. 

In the course of hard photon emission, the value of the 
quantum parameter Xe is shared between the positron 
and the emitted photon [35[, Xe ~ X-y + Xe- If Xe ^ 1 
then both x-y & n d x'e are ^ ess than Xe but are of compa- 
rable value x'e ~ X-y iS Xe- Although propagation of the 
resulting hard photon is not affected by the field, nev- 
ertheless its x-y continues to increase after emission just 
due to rotation of the field. 

In order to understand better what must happen af- 
ter the first hard photon emission, let us come back to 



cvmc 2 // 1 / 4 



e ~ me 2 /// 4 



LOt e 



Huj ' 



1/4 



mc 



x ~ n 



3/2 



(13a) 



(13b) 



(13c) 



Under the condition li ^> 1, as is assumed here, we have 
<C 1 and x ^> 1. The latter inequality approves the 
choice of the asymptotic expressions (|5b[) and (|6b[) . In ad- 
dition, we have the following hierarchy of the time scales 
t a cc *C te, which assures that exactly hard photons with 
X 7 > 1 are typically emitted. 

Within the framework of the doubling chain process 
model, the number of pairs (multiplicity of the cascade) 
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FIG. 4: Left plot: Temporal evolution of the quantum dynamical parameter \ e of the primary electron for three independent 
Monte-Carlo simulations. The thick gray curve corresponds to the analytical solution Eq. (|10b|) for Xe(t) m the absence of 
any QED processes. The three other curves (Runl, Run2, and Run 3) are the results of the three independent Monte-Carlo 
simulations with parameters ao — 2 x 10 4 and Hlj = 1 eV. Right plot: The total number of electrons and positrons iV e + e - vs. 
time for the same independent simulations. 



must grow exponentially, 

iV(*)~e r ', r~i~a^y^. (14) 

In the next section, we are checking the estimations (fl~3|) 
and (fl4| by direct Monte- Carlo simulations. 

IV. DESCRIPTION OF MONTE-CARLO 
APPROACH AND NUMERICAL RESULTS 

In our simulations we are using a Monte-Carlo ap- 
proach for the integration of the cascade equations [see 
the Eqs. (|A"lJ) . ([A2])]. We trace the motion of the elec- 
trons and positrons in between the photon emissions clas- 
sically, whereas for hard photons we exploit the ray trac- 
ing approximation in between their emission and con- 
version into pairs. Even though there exists the exact 
analytical solution ([9]) for equations of motion <JSj> for 
positrons and electrons, we are integrating Eq. ((5} nu- 
merically for each of the particles. This is done in order 
to incorporate the probabilistic events of photon emission 
and pair creation in the routines as described below, as 
well as for the purpose of future generalization to more 
realistic field configurations. 

Our numerical algorithm works as follows. At each 
time step U < t < ti + At we are calculating the momenta 
of all the particles created at preceding time steps by 
Pi+i = Pi + q t E. l+1/2 At, where E t+1/2 = E(U + At/2) 
and qt = +e, — e,0 for positrons, electrons and photons, 
respectively. The event generator determines which of 
the electrons or positrons is going to emit a photon at 
this time step and whether any of the present photons is 
going to produce a pair. 

Let us explain our event generator for photon emission 
in more details (see also Refs. [22|, HU). Starting from 



Pi and Ei = E{tj), we attach the value Xi a t time t, 
using Eq. <j2j) to each electron and positron and compute 
the total probability rate W ra d (see Eq. In order to 
isolate the infrared singularity, we set the lower limit of 
integration to For each electron and positron, we 

assume that it emits a photon between ti and if r < 
W ra dAt, where r (0 < r < 1) is a uniformly distributed 
random number. If the above inequality is fulfilled, then 
the energy e 7 of the emitted photon is obtained as the 
root of the sampling equation 

1 J dW rad (e 7 ) = 
W rad J de 7 7 ' V ' 

where r' (0 < r' < 1) is an independent random number. 
The time step At, which remains fixed in the course of 
computation, must be chosen such that At W~ d , W~ r x 
holds. The direction of propagation of the newly emitted 
photon is parallel to the momentum pi of the parental 
electron or positron, whose momentum after emission we 
find from the conservation law as discussed in Sec. ITT] For 
pair creation, the event generator works similarly, apart 
from the fact that there is no need for the regularization 
parameter e m i n - 

Within the constant crossed field approximation ap- 
plied here we assume that £ 7 ^> mc 2 . However, the pho- 
tons with energies e 7 < mc 2 are not able to create pairs 
in a subcritical field, for which x-y 1 holds. Therefore, 
we can completely neglect emission of soft photons in our 
problem. Based on this reasoning, we currently set the 
lower integration limit £ ml „ to mc 2 for both W rad and 
the sampling Eq. (fT5|) . 

As a benchmark for our code we have simulated the 
development of a cascade initiated by a high-energy 
(eo = 2 x 10 5 mc 2 ) initial electron in a constant homo- 
geneous transverse field with Eq = 0.2Es- Our results 
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FIG. 5: Left plot: The dynamical quantum parameter (xe) for the electrons averaged over the cascade vs. time for the same 
simulations as in Fig. [4] Right plot: Evolution of the mean energy of the electrons, positrons, and photons averaged over the 
cascade in a typical simulation run (ao = 5 x 10 4 and ftui = leV). 



are averaged over 10 3 simulation runs. In this particular 
simulation, the curvature of trajectories of electrons and 
positrons has been neglected, so that the results of our 
simulations can be directly compared with previous sim- 
ulations of cascades produced by high-energy electrons 
in a magnetic field 24'] . Comparison of cascade profiles 
obtained in both simulations is given in Fig. [3] by the 
solid red and dashed green lines, respectively. The figure 
represents the number of pairs with an energy exceed- 
ing 0.1% of the energy of the primary electron versus the 
elapsed time. In our notation the reference characteristic 
radiation time t ra d as adopted in Ref. [HI in our notation 

is t rad = 3.85 x (yn/ax 2 ^ 3 ) x (h/mc 2 ) = 5.64 x W'^, 
where the subscript "in" refers to the initial data for pri- 
mary electron. We see that our results are in reasonable 
agreement with the paper p4| . 

We have also implemented and tested a different event 
generator, which provides significant speed up due to the 
absence of numerical integrations. The idea is to exploit 
some explicit algebraic fits for the energy spectrum ([lj, 
and to exchange the order of testing the occurrence of 
photon emission and of sampling its energy. In this al- 
ternative version of the algorithm, within each time step 
one first samples the possible energy of an emitted pho- 
ton just as an uniformly distributed random quantity, 
e 7 = e e r' in the above notation. After that, photon emis- 
sion is assumed to take place if r < [dW ra d(Sj) /d£j]e e At. 
In this case, the time step must satisfy the condition 
At < [e e dW ra( i(e*)/c?e 7 ] _1 for all appearing electrons 
and positrons, where e* is the photon energy that cor- 
responds to the maximum of the emission spectrum ([lj . 
The same scheme can be applied to the simulation of pair 
photoproduction as well. Note that in this case there is 
no need to introduce the energy cutoff e m i n , although 
this may serve as a useful trick if one wants to restrict 
the number of soft photons that arc traced by the code. 
The test of the modified event generator is included by 



a dashed blue line in Fig. [3] This test demonstrates that 
both versions of the event generator are in fact equiva- 
lent. 

The results of our simulations are collected in the fig- 
ures [21 UK Fig. [2] represents a typical spatial picture 
of the formation and development of a cascade initiated 
by a positron. The electrons and positrons are deflected 
by the field in opposite directions, whereas the directions 
of propagation of photons are distributed randomly, as 
could be expected. For the rest of the paper we assume 
for all simulations in an uniformly rotating field that at 
t = we have a single electron at rest (p e = 0) and 
no photons and positrons. The typical evolution of the 
quantum dynamical parameter Xe of the primary electron 
is illustrated with the left panel of Fig. [4] Before the emis- 
sion of a first photon, the electron is gaining energy and 
its parameter Xe is growing as the square of time in ac- 
cordance with Eq. (|llb[) . After the first photon emission, 
which for our parameters happens typically on the time 
scale t e smaller than lu^ 1 , the curves become stochastic 
and consist of smooth sections with typical growth of Xe 
due to acceleration by the field. These sections are sep- 
arated by sudden breakdowns resulting from recoils due 
to successive photon emissions. Since these recoils are 
random, the three curves in the figure corresponding to 
independent simulation runs deviate at later times. Af- 
ter the transient period which typically lasts for several 
lifetimes t e , the momentum losses due to quantum recoils 
are coming into equilibrium on the average with the trend 
of acceleration by the field. After that, function Xe(t) for 
an individual electron describes a stationary stochastic 
process. 

As it was predicted in Ref. [13], the development of a 
cascade results in exponential growth with time of the 
total numbers of secondary hard photons and electron- 
positron pairs. This is illustrated with the right panel in 
Fig. 2J The plot N e - e + (t) is given by a random stairway, 
with each stair corresponding to creation of a single pair. 
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The successive stairs are well separated initially, when 
the total number of pairs remains small. At later time 
with the number of pairs growing rapidly the stair-like 
structure of the lines in the plot becomes invisible and 
straight lines are obtained. Although these straight lines 
for independent simulation runs are typically different, 
mostly because emission of the first photon starts ran- 
domly from one simulation run to another, nevertheless 
their gradients are varying weakly in different runs and 
can be used to determine the growth rate T in Eq. (fT4l . 
For example, the growth rates extracted from the curves 
1-3 at Fig. Hare V = 4.62, 4.84 and 4.90, respectively. 

We have studied the averages of the quantities Xe, £e 
and over the cascade. For example, temporal evolution 
of the mean value 



10" 



(Xe(t)) = 



1 



N 



N e -(t) 



(t) 



(16) 



where N e - (t) is the instant number of present electrons 
and Xei(t) is the instant value of the quantum dynami- 
cal parameter for the i-th electron, as depicted in the left 
plot of Fig. One can see that at later times the ran- 
dom fluctuations are smoothed out and the quantity (|16p 
stabilizes acquiring some definite constant value which is 
independent of the simulation run. The same behavior 
was observed for (s e ) and (9), which are defined in a 
manner similar to Eq. (|16p . The typical evolution of the 
averaged energy of all the components of the cascade is 
represented in the right plot of Fig. [5] At later times, 
the mean energies of electrons and positrons coincide as 
is expected from symmetry consideration, whereas the 
mean photon energy typically remains smaller. At the 
same time, the energy spectrum of created electrons and 
positrons is wider than the photon spectrum (see Fig.[S]). 
Both features are explained naturally by the fact that in 
our setup the hard photons (\j ^ 1) are quickly con- 
verted into pairs which survive, whereas soft photons 
(x-y ^ 1) ar c stable with respect to pair photoproduction 
and hence are accumulated. In the high energy region, 
all the spectra are likely to show exponential decay. 

One of our main tasks was the investigation of the va- 
lidity of estimations (|T3l) and ((14)) which were suggested 
previously in Ref. [10( and are of crucial importance. In 
particular, Eq. (fT4|) was serving as an argument for the 
estimation of the maximal value of the intensity attain- 
able with focused laser fields. In order to test these es- 
timations we have performed parametric studies of the 
stabilized values of the quantities (xe), (£e), {&) and of 
the increment T. These results are presented in Fig. [7] 
The left plot demonstrates the dependence of the ratios 
of the quantities obtained by simulations to their estima- 
tions ([T3]) on fi for fixed rotation frequency Huj = leV. 
It is clear from the figure that for large values of \i each 
ratio acquires some definite limit of the order of unity. 
According to our results the formulas (|13|) arc valid up 
to some numerical coefficients of the order of unity, which 
vary no more than twice in the whole range /U > 1. The 
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FIG. 6: The energy spectra for different components of the 
cascade at t = 1.2 x lj^ 1 for ao = 5 x 10 4 and huj = 1 eV. 

results of simulations for T are compared with Eq. (Tl"4"|) 
for two different rotation frequencies Huj = 0.66eV and 
Huj = leV on the right panel of Fig. [7] One can see that 
for large values of /i the estimation (fT4|) is justified with 
good accuracy even without any correction factor. For 
fi < 30 formula (|14[) overestimates T but not more than 
by half of an order of magnitude. This may be neverthe- 
less crucial for the estimation of the total cascade yield 
due to its exponential dependence on T. For the par- 
ticular value /j, f=s 10, which was exploited in Ref. [lOj, 
formula (|14|) overestimates T by approximately a factor 
of 1.5. This, however, can be compensated in principle 
by simultaneous underestimation of the escape time in 
Ref. (13 . 

In order to apply the results of our simulations to es- 
timate the cascade yield by a realistic focused laser field, 
we assume that the appearing electrons and positrons 
are pushed away as a whole from the focus by the pon- 
deromotive potential in radial direction with almost the 
speed of light. Assuming the Gaussian profile of the fo- 



cused beam we can write fj,(t) 



Hoe 



where no 



is the value of the parameter \i at the center of the focus 
and wo is the focal radius. Then the total number of 
pairs produced by the cascade can be estimated to the 



ln(iV e+e -)~ / T(ix(t))dt = T(n ) / e 



The remaining integral defines the effective time of escape 
from the focus and equals v / 7r(wo/c), i.e. is y/n w 1.77 
times larger than it was assumed in Ref. [Toj . This cor- 
rection almost totally cancels the overcstimation of T by 
formula (fT4"|) at /i w 10 that we have observed in our sim- 
ulations. Thus, we hope that the quantitative predictions 
in Ref. (loj must remain unaffected by our corrections. 

One can see that we are currently neglecting the com- 
plicating details in our code such as, e.g. the elastic colli- 
sions, the Compton scattering, and the annihilation pro- 
cesses. Such phenomena must become important only at 
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FIG. 7: Left plot: Parametric studies of the mean energy e e , the mean dynamical quantum parameter Xe> and the mean angle 
6 e between the momentum and the field for electrons and positrons. The ratios of the simulation results to the approximations 
()13bl) . (|13c|) are plotted vs. the parameter fi for fuo = leV. Right plot: Parametric study of the increment F as a function of 
the dimensionless field strength \x for two rotation frequencies ftw = 1 eV and fku — 0.66 eV. The approximation (|14[) is shown 
by the dashed lines. 



later times, when the plasma is dense enough. Though 
successive collisions and annihilations of the electron and 
jositron from the same created pair may be important 
30 1 ] , we are currently ignoring these effects for simplicity 
3q . We ignore the higher-order processes, such as two- 
photon creation and the trident processes as well (see the 
remark at the end of Sec.[II|. All these assumptions are 
natural and commonly accepted in the present cascade 
theory [23|, even though they may be revised in future 
studies. 

Due to limitations of computer power, we cur- 
rently stop our simulations after the creation of around 
N e - e + < 10 4 pairs. This was shown to be enough to 
estimate the growth rate T, as well as to average the 
characteristics of a cascade over the ensemble of pairs 
with reasonable accuracy. In the time interval of simula- 
tion these pairs occupy a volume of the order d 3 , where 
d ~ c/ui ~ 1/im. This corresponds to the pair density 
n e + e - ~ 10 16 cm -3 . Typical values of the 7-factor for 
electrons and positrons arc 7 ~ 10 3 ™10 4 , see Fig. [6] This 
corresponds to their energies e e = jmc 2 ~ 0.5 5GcV. 
Assuming the temperature T ~ e e /k ~ 10 13 K, we can es- 
timate the Debye screening radius rrj ~ \JkT '/ 'e 2 n e + e - ~ 
lcm 3> d. The relativistic plasma frequency VL pe ~ 
c/rjj ~ 10 sec remains about five orders of magni- 
tude smaller than the optical frequency. For these rea- 
sons we have completely neglected Coulomb interaction 
between electrons and positrons and all the accompany- 
ing collective plasma effects in the present simulations. 
However, the density of pairs is growing exponentially, 
n e + e - (t) cx e rt , and hence ro oc e~ Tt / 2 and Q pe oc e rt / 2 . 
After a relatively short period of time < 2ir/ui, when 
the number of pairs becomes macroscopic (~ 10 11 ), the 
quantities tjj and Q pe attain the values d and w, respec- 
tively, and the collective plasma effects may come into 
play. Within our approximation of a homogeneous field, 



the total number of created pairs would be restricted by 
the screening of the external field by the self-field of aris- 
ing plasma. 

Let us note, that despite some doubts in the literature 
[25l | . the radiation friction is taken into account properly 
in our version of the algorithm by the recoils happen- 
ing at the times of photon emission (see, e.g., Ref. (3l| 
and the Appendix lAl in our paper). Hence, there is no 
need to include an additional radiation friction force in 
the equations of motion (|5J) for electrons and positrons, 
otherwise this would cause double counting. Moreover, 
our approach transfers the concept of classical radiation 
friction into the quantum domain in a correct fashion. It 
can be asked how the classical continuously acting radi- 
ation friction can be recovered from the sudden jumps 
of momentum similar to those in Fig. |H In fact this 
happens on the average with respect to the ensemble of 
Monte-Carlo realizations, since the moments of succes- 
sive photon emissions arc distributed randomly. At later 
times, when the number of created pairs becomes large, 
the cascade forms a representative ensemble itself and 
there is in principle no need for taking the average over 
independent realizations. 

V. SUMMARY AND DISCUSSION 

In this paper we have presented the first results of nu- 
merical simulations of the formation and development of 
electron-positron-photon cascades by initially slow elec- 
trons in a uniformly rotating homogeneous electric field. 
In such a situation the cascades reveal a principally new 
feature, i.e. the restoration of the energy and the dy- 
namical quantum parameter due to the acceleration of 
electrons and positrons by the field. This feature may 
be of crucial importance for the whole physics of laser- 
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matter interactions in the strong field domain, as it was 
demonstrated in Ref. [l(|. We have explicitly identified 
this restoration mechanism in the course of our simula- 
tions. Also, our simulations clearly confirm the qualita- 
tive analysis of Ref. [Io[ , including the basic scalings (fT3|) 
and the estimation (Q3J) for the cascade yield. So, they 
can be used to fix the remaining numerical correction fac- 
tors in Eqs. f) 13[) and (HH), which turn out to be of the 
order of unity. 

The numerical approach that we adopt is based on 
Monte-Carlo simulations of the cascade equations. We 
have shown explicitly that contrary to some recent 
doubts in literature [25] such an approach incorporates 
radiation friction acting on individual electrons and 
positrons and, moreover, is doing this in a manner which 
is consistent with intense field QED. 

The codes designed for our task can be readily adopted 
for simulating cascades in the laser fields with more re- 
alistic configurations, such as tightly focused Gaussian 
beams and pulses. This is required in order to make 
more definite predictions on the impact of cascade pro- 
duction for possible future experiments, as well as for 
further corrections of the maximal value of intensity that 
can be attained with optical lasers floj |. Simulation of 
cascades in a focused laser field will be presented in a 
separate publication. However, let us make several brief 
comments about cascades in focused laser fields, possible 
experimental scenarios, and some yet unresolved techni- 
cal problems that may require further studies. 

The restoration mechanism arises due to the curva- 
ture of the trajectories of the charged particles across the 
field and may be sensitive to its polarization. Although 
we expect that this mechanism must work for generic 
field configurations (e.g., for generic tightly focused laser 
fields), there may exist several particular configurations 
for which the restoration mechanism does not work. For 
example, in an arbitrary constant electromagnetic field or 
a circularly polarized propagating plane electromagnetic 
wave the dynamical quantum parameter Xe is conserved 
exactly in the course of motion. In the case of a generic 
propagating plane wave the amplitude of oscillations of 
the parameter \e for an initially slow electron does not 
exceed Eq/E$, i.e., always remains smaller than unity. 
Another example is a linearly polarized oscillating elec- 
tric field [25|, since in this case the initially slow parti- 
cles are accelerated strictly along the field and hence the 
growth of the transverse component of the field is absent. 
In some intermediate cases, e.g. for elliptical polarization 
or a weakly focused Gaussian beam, restoration of Xe 
must exist but may be less effective than in the case of cir- 
cular polarization. However, in all the cases at least the 
usual cascades would be caused by external high-energy 
electrons or hard photons passing through the high-field 
region transverse to the field. In this case, the cascade 
yield remains microscopic and would be determined by 
both the initial energy of an external energetic particle 
and the laser field strength. 



In order to initiate a cascade in a tightly focused laser 
field, it is required to inject a primary particle into the 
center of a focal region. This task may be not triv- 
ial because the focal region is surrounded by a pon- 
deromotive potent ial wall of the characteristic height 
Uq ~ mc 2 yl + a 2 , w mc 2 ao- The external high-energy 
electrons will be most likely deflected rather than pen- 
etrate inside. The most direct and elegant scenario is 
based on the exploration of pairs that are created spon- 
taneously from vacuum by the laser field itself [l(| , since 
they are appearing exactly at the center of the focus as 
required. However, this possibility implies high intensi- 
ties > 10 26 W/cm 2 . Another possible resolution would 
be the initiation of cascades by energetic 7-quanta. In 
our opinion, the final conclusion whether or not cascades 
with macroscopic yield can arise in generic real exper- 
iments exploring laser-matter interaction at intensities 
lower than ~ 10 26 W/cm requires further studies. We 
note that for cascades that arise in the course of the in- 
teraction of high-intensity laser radiation with material 
targets it may be necessary to take the impact of ordinary 
cascades in matter into account as well [331 ] - 

If the cascade yield attains macroscopic values 
(N e - e+ ~ 10 11 ), the self-field of the electron-positron 
plasma becomes comparable to the guiding field. In this 
regime screening of the external field and its absorption 
by the electron-positron plasma self-field will restrict fur- 
ther pair production. Such a regime can be simulated 
by combining our codes with the Particle-in-Cell (PIC) 
method [Hj]. We hope to address this problem in one of 
our next publications. 
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Appendix A: Cascade equations and radiation 
reaction 

The cascade equations for a uniformly rotating homo- 
geneous electric field 
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df±(Pe,t) ± eE , t s df±(p e ,t) 



at 



Ope 



d/ 7 (p 7 ,Q 
dt 



= J Wrad{Pe+P~f ^ P 1 )f±{Pe+P~ n t) d 3 p~ f - f±(p e ,t) J W rad (Pe ~> P 7 ) d 3 p 7 

+ J w CT (pj -*-p e )/ 7 (p 7 ,i)d 3 p 7 ,(Al) 

J W rad (p e -> P 7 )[/+(p e ,t) + f-(Pe,t)]d 3 p e - / 7 (p 7 ,t) J W cr {p 1 -> p e )d 3 p e .(A2) 

I 



differ from the standard equations of EAS [12|,|13| only by 
the addition of the second term to the LHS of Eq. (|Al|) , 
which takes electron and positron acceleration into ac- 
count. Here, f± and / 7 are the distribution functions 
for positrons, electrons and photons, respectively. In our 
approximation of photon and pair emission in strictly for- 
ward direction the differential probability rates are of the 
form 



Wrad(Pe ~> P 7 ) = J d 



dW raa 



£ 7 — Ae e 



Wcr(Py -> Pe) = / 



cfe e 



5(p e - Ap 7 )(A3) 



so that the integrals standing on the RHS of Eqs. (|A1[) . 
(|A2p are essentially one-fold. Nevertheless, in the main 
case of interest, the direction of the field E(t) varies in 
time, so that the problem does not reduce to ID. Note 
that the scalings (fT5|) can be obtained from (|A1[) and 
(|A2p in the limit of large x y i a dimensional analysis. 

It is interesting to note that in the view of the approx- 
imation (|A"3)) it follows from the Eqs. (|A"Tj) . (|A"2|) that 



A i / ^ 



+ /_)d 3 p e + / P 1 f 1 d 3 p 1 



dt 



eE(t) / (/+ - f-)d 3 p e , (A4) 



(A5) 



e e (/+ + f-)dp e + J e 1 f 1 d 3 p 1 
= eE(t) f El(f + -f_)d 3 Pe 



These identities ensure that in our approximation the 
momentum and the energy are extracted from the field 
only during the acceleration of electrons and positrons, 
i.e., both energy and momentum of the clcctron-positron- 
photon plasma arc conserved during photon emission and 
photons to pairs conversion. 

The first two terms on the RHS of the Eq. (jAl)l de- 
scribe the influence of photon emission on electrons and 
positrons. Let us now demonstrate that classical radi- 
ation reaction is taken into account properly by these 
terms. 

For this aim let us assume in what follows that the 
electrons are slow in the sense that distributions f± are 



restricted to such momenta for which Xe± "C 1 (though 
we continue to assume that they are ultrarelativistic). 
In this case their motion can be described in completely 
classical terms. Accordingly, let us skip the third term 
on the RHS of the Eq. (|A1|) . which is responsible for 
pair production. Once this is done the total numbers of 
positrons and electrons N± = J f±d 3 p e are conserved. 

The relation between the variables % 7 and x in Eq. ([1]) 
can be expressed in the alternative form 



X 7 



x 3,2 xl 
1 + x 3 / 2 x 



(A6) 



It is clear from this formula by taking into account that 
the spectrum ((T|) of emitted photons is effectively con- 
centrated in the range x < 1, that 



X-y 



C 3/2 Xe < Xl < Xe 



(A7) 



As a consequence, p 7 <C p e in the remaining integrals on 
the RHS of Eq. ((ATI) . Thus, the expansion 



Wrad(Pe + P 7 



P 7 )/±(Pe) 



P 7 )/±(Pe + P 7 ) - W rad (p e 

d 

~ Pi ■ a \ W rad{Pe ~> P 7 )/±(Pe)], 



dp e 



(A8) 



is valid. 

Consider the average momenta P±{t) = 
(l/N ± )Jp e f ± (p e ,t)d 3 p e . By multiplying Eq. (JAT|, 
truncated as described above, by p e and integrating it 
over p e by parts we have 



P±(t) = ±eE(t) + (R) ± (t), (A9) 



where 



R(p e ) = - J P 7 «W(Pe P 7 ) d 3 p~ f (A10) 

is the mean rate of momentum losses of electrons and 
positrons due to photon emission, and 

(R)±(t) = jLJ R(p e )f ± (p e ,t) 

are its averages over the momentum distributions of 
positrons and electrons, respectively. 
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Taking into account that p 7 « x 3 / 2 \ePe and e 7 = 
(e e /Xe)x 7 and using Eqs. (|A3|) and |T|), we arc coming 
to 

Xe 

H(p e ) = -e ePe J x^^dx,. (All) 



From this point, let us pass to the integration over the 
variable x. In the approximation \e <C 1 we have x 7 = 
x 3//2 Xg. Contribution to the integral in Eq. (|Alip comes 
from the range x ~ 1. We can neglect the term x-yV% 
in the brackets in the expression ([I]) and, in addition, 
replace the upper limit of integration over x by infinity. 
After these manipulations, we have 

o 9 4 2 

6am z c*xiPe T 

H = o fc J ' 

2 ft £ e 

oo f oo "\ 

J = J \x 2 J Ai(f) d£ + 2a;Ai , (.T) I da. (A12) 

The remaining integral J by integration by parts and ex- 
ploiting the Airy equation Ai (x)-xAi(x) = is reduced 
to 

oo 

2 f 4 
J = — - / x 3 Ai(a;) = --. 

3 J 9 

o 



Thus, in a view of Eq. (J2j) , we finally have 



R 2 am 2 c 4 Xe Pe 2 e 4 f I7 2 cp e 
3 ft £ e 3 m 2 c 4 e e 



This is exactly the main contribution to the Landau- 
Lifshitz (LL) force [26| for ultrarelativistic electrons. 
Other contributions to LL equation do not appear in our 
derivation only because we are essentially using the con- 
stant crossed field approximation. 

In the quantum case %e 3> 1 the expansion (|A8j) is not 
valid. Thus, radiation friction in the quantum regime 
cannot be described by the concept of classical force in 
principle, as it was attempted to do e.g. in Q. In ad- 
dition to the advection term in the transport equation, 
which could be ascribed to the radiation reaction force as 
above, spreading of the distribution functions in momen- 
tum space becomes important as well. This spreading is 
associated with quantum fluctuations and is observable, 
e.g., as quantum excitation of synchrotron and betatron 
oscillations [HJ. 
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